Finite temperature coherence of the ideal Bose gas in an optical lattice 
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Abstract 

In current experiments with cold quantum gases in periodic potentials, interference fringe contrast is typically the easiest 
signal in which to look for effects of non-trivial many-body dynamics. In order better to calibrate such measurements, we 
analyse the background effect of thermal decoherence as it occurs in the absence of dynamical interparticle interactions. We 
study the effect of optical lattice potentials, as experimentally applied, on the condensed fraction of a non-interacting Bose 
gas in local thermal equilibrium at finite temperatures. We show that the experimentally observed decrease of the condensate 
fraction in the presence of the lattice can be attributed, up to a threshold lattice height, purely to ideal gas thermodynamics; 
conversely we confirm that sharper decreases in first-order coherence observed in stronger lattices are indeed attributable to 
many-body physics. Our results also suggest that the fringe visibility 'kinks' observed in F.Gerbier et al., Phys. Rev. Lett. 
95, 050404 (2005) may be explained in terms of the competition between increasing lattice strength and increasing mean gas 
density, as the gaussian profile of the red-detuned lattice lasers also increases the effective strength of the harmonic trap. 



I. INTRODUCTION 



The first investigations of dilute Bose-Einstein conden- 
sates (BECs) in periodic potentials [l|, 0, Q launched 
a research field which continues to expand steadily 
[J, Q. The theoretical interest of the subject Hes in 
the wide range of analogous phenomena for which it 
can provide an idealized model: spin-spin interactions 
in crystalsjg, 0| and general Bose-Hubbard Hamiltonian 
dynamics|8|, effective magnetic fields with non-abelian 
gauge potentials IqI, [ig|. and even gravitational phenom- 
ena [ll|. Experiments exploiting the high tunability of 
optical lattices have observed the Mott insulator phase 
transition (l2 . [l3J . the Berezinskii-Kosterlitz-Thouless 
transition [l^, strong quantum depletion [l5|, extended 
coherency with number-squeezed states [l6|, second-order 
atom tunneling [l7[, resonantly enhanced tunneling 18], 
and repulsively bound atom pairs [19[. The combination 
of tunable optical lattices and ultracold temperatures 
provides enough control over atomic motion to investi- 
gate numerous questions of fundamental interest and im- 
portance. 

A workhorse experimental technique for such investi- 
gations is the measurement of matter wave decoherence. 
Measurements of collective excitation spectra can pro- 
vide independent information about many-body dynam- 
ics [12i | ; but such measurements are in general somewhat 
more involved. The coherent population fraction can be 
inferred fairly straightforwardly from measurements of 
fringe contrast in matter wave interference. A drop in 
coherent population fraction can signal the formation of 
interparticle correlations due to many-body interactions 
[201 ■ It can also be caused, however, by simple heating 
of the sample. Moreover, ramping up an optical lattice 
alters the density of states in a gas, and so can change its 
equilibrium properties even if heating as such can be sup- 
pressed. To use coherent population fraction as a signal 
of many-body physics, therefore, one must understand 
these background effects well enough to subtract them 
accurately. 

The question of temperature and critical tempera- 



ture change during the process of loading an optical lat- 
tice with atoms has already received substantial study 
[H m, [13, m, ii, [13, il III . Unfortunately the prob- 
lem turns out to be quite complicated in practice. While 
it is straightforward to show that ramping up a uni- 
form lattice extremely slowly will adiabatically cool a 
non-interacting ultracold gas, interactions and a time- 
dependent trapping potential are unavoidable experimen- 
tal complications. And although experiments are able to 
allow quite long times for the gas to equilibrate after the 
lattice has been fully turned on, the available experimen- 
tal lifetimes set limits on how slowly the lattice potential 
can be ramped up, in order to reach full lattice strength 
before the gas escapes the trap. Under present conditions 
it therefore seems difficult to ensure that the gas really re- 
mains in equilibrium during the entire experiment. Some 
loss of condensate fraction does occur when the lattice 
strength is raised and lowered in a cycle, and although 
it is small enough to allow ample observations of Bose- 
Einstein condensates in optical lattices, it does indicate 
non-equilibrium entropy generation during ramping. 

The theory of heating during loading of a condensate 
into an optical lattice is thus a subject of active research. 
It is an important task for cold atom theory, because di- 
rect experimental measurement of very low temperatures 
in highly condensed Bose gases remains technically diffi- 
cult. (The standard fitting of the thermal fraction's mo- 
mentum distribution to the Maxwell-Boltzmann formula 
becomes problematic when the thermal fraction is small.) 
The purpose of this paper is to study the effects of lattice 
potentials, as experimentally applied with focused laser 
beams and with harmonic trapping potentials superim- 
posed, on the density of states of the ideal Bose gas. Our 
main contribution is to identify a basic competition be- 
tween density of states and spatial density, which occurs 
because of the particular way optical lattices are imple- 
mented in current experiments, and which considerably 
moderates the effect of weaker lattices on gas coherence. 

As a first exercise in identifying the quantitative conse- 
quences of these effects in experiments, we model several 
recent experiments under the simplifying assumptions 
that the gas does remain in global thermal equilibrium 



at all times, and moreover that it maintains a constant 
temperature. We determine this temperature by fitting 
to experimental data from before the lattice potentials 
are ramped up. This scenario would be expected in the 
case where the trapped gas was maintained in equilibrium 
through scattering with a reservoir of thermal gas, which 
was itself too slightly perturbed by the time-dependent 
lattice to undergo any temperature change. For shallow 
lattices our isothermal model turns out to fit experimen- 
tal data somewhat better than the alternative scenario 
of isentropic evolution optical lattices, which would be 
expected if the lattice-trapped gas could be regarded as 
isolated, and if the lattice were ramped up slowly enough 
for equilibrium still to be maintained. Our quantitative 
results thus provide some somewhat surprising input for 
the further non-equilibrium studies which will ultimately 
be required in this problem. Moreover, they serve to con- 
firm the dynamical competition effect we identify, which 
must play an important role in determining how the gas 
coherence depends on the lattice potential, regardless of 
the precise statistical ensemble which should describe the 
gas. 

The paper is organized as follows. We first indicate the 
parameter regime in which an ideal Bose gas in thermal 
equilibrium is expected to be a good model for exper- 
imental systems. We then describe the theory we will 
use for a Bose-condensed ideal gas in an optical lattice 
at finite temperature. Along with this full version of 
our ideal gas theory, whose predictions we will compute 
numerically, we introduce some further approximations 
that will allow analytical results, and help interpret our 
numerical curves. 

Wc then present our numerical calculations and com- 
pare them with several sets of experimental data, includ- 
ing experimental configurations in which the gas dynam- 
ics is effectively three-, two-, or one-dimensional. In par- 
ticular we consider the fringe visibility "kinks" seen in 
recent experiments, where the downward sloping curve 
of first order matter coherence versus increasing lattice 
strength has a shoulder. We show and explain a very 
similar effect in the ideal gas coherent fraction, which 
occurs at closely comparable lattice strengths. 

In our final section we summarize our results and pro- 
vide an outlook sketch. 



II. THERMAL EQUILIBRIUM 

We assume interaction between particles sufficient to 
maintain equilibrium by Boltzmann scattering. In a 
truly coUisionless regime, an adiabatic change of poten- 
tial would simply preserve the population in each energy 
level, and so would not change the number of condensed 
atoms. In considering an ideal gas in equilibrium, how- 
ever, we naturally mean that we consider time scales long 
enough for interactions in the form of collisions to gener- 
ate the time-averaged behavior that is represented by a 
grand canonical ensemble; but we do not admit interac- 



tions strong enough to alter the instantaneous probability 
distributions in that ensemble, which therefore depends 
only on the non-interacting Bose gas Hamiltonian. 

Cold dilute gases in experiments can certainly be made 
very weakly interacting, so the experimental relevance of 
our model depends on the relation between the experi- 
mental time scale tx and the timescale tc for coUisional 
equilibration. Quantum kinetic theories for ultracold 
atoms support the identification of the latter timescale 
with Boltzmann's scattering time: 
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where p is the gas density, T is the gas temperature, a is 
the s-wave scattering length for low momentum transfer 
collisions, and m is the particle mass. This gives an equi- 
libration timescale on the order of tc ~ AOms for ^^i?6 



with densities of p w 10 



Since lattice ramp time 



scales in experiments are generally much longer than this, 
our assumption of equilibrium because tc ^ tx should 
be adequate for experimental comparison. 

The validity of our ideal gas approximation then de- 
pends on the ratio between tc and the relevant time scale 
T]j for single-particle dynamics. What then is this td? 
A global definition in terms of the rate of change of sin- 
gle particle energy eigenstates would be appropriate for 
experiments in the coUisionless kinetic regime, studying 
the first onset of interaction effects. For current exper- 
iments approaching phenomena such as quantum phase 
transitions due to onsite interactions, a local criterion is 
more relevant. We therefore consider td to be the time 
scale for quantum tunneling between two adjacent optical 
lattice sites. We can invoke the energy-time uncertainty 
relation to estimate this tunneling time as td — h/SE 
for SE the width of the lowest energy band. At opti- 
cal lattice depths of order 40Er {Er = h^k'^ /2m being 
the 'recoil energy') this implies T]j « lOOms. For the 
strongest experimental lattices we thus have td ^ tc, 
so that particles typically collide before leaving one lat- 
tice site. The ideal gas approximation therefore breaks 
down, and interparticle interactions begin to have non- 
trivial effects beyond merely enforcing equilibrium over 
long time scales. Accurately assessing such effects is one 
of the major challenges in current many body theory; our 
goal in this paper is simply to identify the onset of this 
regime in current experiments. 

For shallower optical lattices, however, the tunneling 
time rapidly decreases: td ~ 1ms for a lattice depth of 
12Er. With Td <S^ Tc <S^ Tx , particles will typically tun- 
nel through many lattice sites before scattering, so that 
the resulting equilibrium distribution of single-particle 
energies will include the band structure of the lattice, but 
the effects of harmonic trapping over length scales much 
longer than the lattice spacing can be computed in local 
density approximation. This theory should apply up to 
the onset of dynamically (as opposed to kinetically) sig- 
nificant interactions at high lattice potentials. As we will 
see, it does indeed fit experiments well up to threshold 



lattice strengths. This agreement is not trivial, however, 
but depends upon the cancellation of two basic features 
of the single-particle dynamics in typical experimental 
potentials. The interaction onsets identified in this pa- 
per therefore provide useful information for interpreting 
experiments on ultracold atoms in optical lattices. 



III. MODEL POTENTIAL FOR OPTICAL LAT- 
TICE PLUS HARMONIC TRAP 

We consider an ideal Bose gas in a harmonic trap to 
which has been added a periodic potential created by 
counterpropagating laser beam pairs in a;, y and z di- 
rections. It turns out to be very important to represent 
the experimental fact that the focused red-detuned laser 
beams which generate the lattice potential have gaussian 
radial envelope profiles [30J: 
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Here a;o,i is the angular frequency of the harmonic trap- 
ping potential in the zth direction, k is the wavenumber 
of the ID lattice potential formed by the AC Stark effect 
in the standing waves of the ith direction's laser pair, and 
Wi is these beams' 'waist' (radius at which laser intensity 
has fallen by a factor 1/e^). We write s\r?{kri) — 1 in- 
stead of — cos^(fcri) because we will treat differently the 
smooth gaussian envelope of the lattice minima and the 
gaussian modulation of the lattice barrier height. 

The importance of the finite waist in the gaussian beam 
profile is as follows. All of these experiments use red- 
detuned light, so that the laser intensity maxima provide 
the lattice potential minima. This means that the AC 
Stark potential felt by the atoms is equivalent to a broad 
gaussian well, plus a repulsive lattice whose height dimin- 
ishes with radius (according to the gaussian envelope). 
The spatial extent of the gas cloud is limited by the com- 
bined slowly- varying potential which is the superposition 
of the confining harmonic trap, plus the broad gaussian 
envelope of the attractive laser intensity maxima. Since 
the beam waist is in our cases much wider than the ex- 
tent of the gas cloud in the trap, we can Taylor-expand 
the Gaussian profile of the beam, so that the envelope 
is approximately quadratic over the extent of the cloud. 
We must therefore consider that strengthening the lat- 
tice potential simultaneously strengthens the harmonic 
trapping confinement of the gas. 

As well as the lattice potential minima becoming less 
deep towards the edges of the cloud, the height of the lat- 
tice barriers also decreases with radius; and this changes 
the local density of states for the atoms. For the experi- 
mental cases considered in this paper, however, the vari- 
ation in the local width of the first energy band over the 
extent of the gas cloud is never greater than about 5%. 
This difference becomes even smaller for the stronger lat- 
tices that induce significant decoherence. We will there- 
fore simplify our calculations considerably, at small cost 



in accuracy, by neglecting the effectively small local vari- 
ation in lattice barrier height. The reason we can neglect 
the beam profile for the lattice barriers, but must include 
it in the attractive potential envelope, is that the relevant 
energy scale for the confining effect of the lattice mini- 
mum envelope is the gas chemical potential /i, while the 
energy scale relevant to the effect of lattice height varia- 
tion on bandwidth is at least Er, which is much higher 
than /i in the cases we consider. We will therefore repre- 
sent the potential induced by the red-detuned laser light 
as a superposition of a smooth harmonic well and a uni- 
form repulsive lattice: 
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for i,j, I in cyclic order. 

As a further technical simplification with negligible im- 
pact on the accuracy of our results, we will then replace 
the actual sinusoidal lattice with a parabolic spline po- 
tential, which is piecewise quadratic, and in particular is 
a periodic succession of upward and downward parabo- 
las ('biparabolic'). As Fig.l shows, this approximation 
is very accurate, but it will allow us to compute band 
structures more easily by taking advantage of analytical 
properties of confluent hypergeometric functions, instead 
of dealing with Mathieu functions. We will therefore fi- 
nally write 
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having defined the periodic function 
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where the integer n is chosen for given ^ such that (n — 

l/2)TT^^i^{n+l/2)TT. 

For the uniform biparabolic lattice, {i.e. with uji set to 
zero in V above), the single particle Schrodinger equation 
factorizes and the dispersion relation between energy and 
quasimomentum p is given by a some of one-dimensional 



energies. 
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We show in the Appendix that the Ei{pi) can be found 
implicitly from dimensionless equations of the form 
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where G is a particular combination of confluent hyperge- 
ometric functions [3l| (see the Appendix). These equa- 
tions can be solved numerically to produce an explicit 
numerical interpolating function for E(j)) in the uniform 
lattice, for any given lattice strengths Vi. 




tain 



FIG. 1: Biparabolic periodic potential (solid line) used in this 
paper, in place of the actual sinusoidal lattice (dotted line). 



IV. THERMAL EQUILIBRIUM IN LOCAL DEN- 
SITY APPROXIMATION 

We assume that the state of the (dynamically) non- 
interacting gas may be represented with a grand canon- 
ical ensemble, and this determines the number of non- 
condensed atoms Nnc as a function of chemical potential 
H and temperature T. We also take note of the fact that 
the gas density, locally averaged over a range of lattice 
cells, varies very slowly on the lattice scale. This allows 
us to use a local density approximation, as follows. We 
will consider 
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for fixed f to be the 'local chemical potential', and com- 
pute the grand canonical density of non-condensed atoms 
p{r) in the uniform lattice at temperature T and chemi- 
cal potential /i(r). We will then integrate over r to obtain 
the total Nnc- In the next Section we will then compare 
the theoretical A^atc to the experimentally measured to- 
tal particle number N, for various recent experiments. 
The difference Nc = N — Njyc is inferred as the number 
of ideally Bose-condcnsed atoms. The coherency fraction 
Nc /N will then be interpreted as the coherency fraction 
that would be expected in those experiments if only ideal 
gas physics were involved. 
We therefore write 
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where 1/(3 = fc^T. We perform the p integrals using the 
so-called extended-zone scheme, integrating over quasi- 
momenta within the successive energy bands. In most 
cases we will be at low enough temperatures to restrict 
our integration to the first energy band only. 

To now integrate over r, we expand the Bose-Einstein 
denominator in a rapidly converging Taylor series to ob- 



Nnc = j Srp 
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and, since we are always considering samples in which 
some condensate is present, 

3 

M^^i?.(0). 

i=l 

We can perform the p integrals in the last line of Eqn. (fTU)) 
numerically, for each n, using the E{p) determined by 
numerically solving ^ . The sum converges rapidly, and 
so we obtain a numerical figure for N^c for any given 
temperature T < Tc and set of lattice strengths Vi. 
Eqn. pU]) , evaluated in this way, is the main result of this 
paper. In combination with experimental measurements 
of total particle number, and temperature estimates, it 
will provide the 'theoretical' curves plotted against ex- 
perimental data in the next Section. 



A. Further approximations 

Since our main result (jlOp is evaluated numerically, 
it is useful to compare it with simpler approximations 
whose physical content is more transparent. We will now 
present several of these, of which two will also be dis- 
played, as dashed curves, in some of the next Section's 
plots. 

As the periodic potential becomes stronger, the right- 
hand-side (RHS) of Eq.([7|) approaches a linear function 
of Ei/ Er (within each energy band). As an approxima- 
tion that should become valid for high lattice barriers, 
therefore, we take the RHS to be an exactly linear func- 
tion and obtain the following dispersion relation for the 
first energy band: 
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Here £;""" and E'^^'^ are the first energy band's boundary 
values. Although we determine these band edges from 
numerical calculations for our model lattice potential, 
this approximation for Ei is of the same form obtained 



in the well known tight binding approximation. We will 
therefore refer to this as the tight-binding approximation 
for E{p). Using it, the p integral in (fTO| can be performed 
analytically over the first band, yielding 
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where /q is the modified Bessel function. Eqn. (1121) . ^^~ 
ferred to as 'the tight-binding approximation' for N^c, 
will be used to generate the dashed curves in several Fig- 
ures of the next Section. 

For readers not intuitively familiar with modified 
Bessel functions, of course, this approximated expression 
is no more transparent than the more accurate numer- 
ical results of Eqn. (fTO|) . To indicate the dependence 
on the energy bandwidth more explicitly, however, we 
can simplify Eqn. (|12p further in the strong lattice limit 
Ai <^ ksT^ where the bandwidth shrinks far below the 
temperature, but the temperature in turn remains well 
below the interband gap. In this limit we approximate 
the modified Bessel function as /o(0) = 1 and obtain 
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where Lis (C) = Sfe°=i ^'^/^* is the polylogarithm func- 
tion. 

This simple formula shows well the qualitative impact 
of the lattice parameters on the noncondensed particle 
number. Strengthening the periodic potential compresses 
the energy bands, raising the density of states and thus 
increasing the noncondensed fraction. Increasing the har- 
monic trapping strength instead favors condensation, and 
hence tends to remove non-condensed atoms. Since the 
beam focus means that strengthening the lattice also 
strengthens the effective harmonic trap, these two effects 
compete as the lattice strength is ramped up. 

This simple competition can also be seen, still using the 
tight-binding dispersion relation (jlip . but in the opposite 
temperture regime A^ >> /c^T. Using the asymptotic 
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^NC being the number of non-condensed bosons in a 
harmonic trap of frequencies oji at temperature T [32l |. 
(C(z) is the Riemann zeta function, and C(3) = 1.202). 
Neither of the two simplified formulas (fTS]) and (fT4ll will 
be compared with data in this paper, or for that matter 
referred to again at all; but they serve to illustrate quali- 
tatively the competition between trap and band compres- 
sion that is exhibited in the tight-binding approximation 
P^ and in the full theory of ([TI 



In the regime of temperatures much below the band- 
width, we expect most atoms to explore only the 
parabolic region of the dispersion relation, near its min- 
imum. This immediately suggests, as an alternative to 
tight binding, the parabolic or effective mass approxima- 
tion, in which the motion of particles in periodic struc- 
tures is described as motion of free particle but with an 
effective mass (which is in general anisotropic): 
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By Taylor expanding and keeping the first two terms on 
both sides of our dispersion relation Eq. ([7|) , we obtain 

^ 2m d ^, ,,,„,, 
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Now approximating Ei — > m*pf /2 and integrating over 
the infinite range of p, we find 
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Referred to as 'the effective mass approximation', this 
simplified formula will provide a second dashed curve for 
comparison in some of the Figures of the next Section. 

Thus for low temperatures as well as for higher ones, 
and whether we assume tight binding or quasi-free mo- 
tion, we see that the non-condensed particle number has 
on one hand a tendency to increase (at the expense of the 
condensate) as the lattice barrier height rises, because 
this raises the effective mass and compresses the energy 
bands. On the other hand the non-condensed number 
also tends to shrink, so that the condensate grows, when 
the harmonic confining potential is strengthened. And 
since experiments with focused red-detuned laser lattices 
make the barrier height and trapping potential rise and 
fall together, though not in direct proportion, a non- 
trivial competition is to be expected. 

For temperatures within the bandwidth A^, the non- 
condensed particle number clearly depends on the details 
of the dispersion relation, and no simple approximation 
to Eqn. (jlOp will be very accurate. Our full numerical 
results for Eqn. (8), however, will still represent essen- 
tially this same competition between barrier height and 
containment. 

Before comparing our results with several recent exper- 
iments, we pause to show that the competition between 
confinement and the barrier lattice is indeed important. 
The effectively increased harmonic confinement provided 
by a stronger red-detuned focused lattice can compensate 
to a considerable degree for the increased barrier height, 
and maintain a much higher condensate fraction than 
would survive if the total harmonic confinement were 
truly held constant. If we have a gas of iV = 1.5 * 10^ 
atoms in a harmonic trap with ujx^Oi^y,Q:'-^z,o — 2tt x 
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FIG. 2: First band energy-quasimomentum relation for the bi- 
parabolic lattice (solid line) , in one dimension, with Vi = AEr . 
The dashed and dotdashed lines give the same dependance in 
the tight binding and effective mass approximations, respec- 
tively. 
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FIG. 3: Condensed fraction as a function of optical lattice 
height with (dashed) and without (solid) inclusion of the 
Gaussian profile of the beams. 



{20,20,20} Hz, and T/T° = 0.3(157ii^) is assumed to 
remain constant, then the condensate vanishes entirely at 
an optical lattice height of 14£'r. Including the increased 
effective confinement, due to a gaussian beam profile with 
1/e^ radii of 120 /zm, changes the behavior completely: 
the condensate fraction fc = {N — Nnc')/N changes by 
less than ten percent over this range of lattice heights. 
See Figure 3. 



V. COMPARISON WITH OBSERVATIONS 

Figures. 4 through 7 show the results of our full numer- 
ical theory described above(isothermal) (Eq. (fTO|) ). plot- 
ted as a solid curve (red online), as well as of the two 
semi-analytic approximations (Eq. p^ for tight binding 
approximation and Eq. (|15[) for effective mass approxi- 
mation) and of the isentropic evolution, when their in- 
clusion is instructive. Superimposed for comparison are 
data points directly from the experimental publications. 



We take the experimentally measured total numbers N 
of trapped particles, together with our theoretical values 
for non-condensed particle numbers A'jvd to compute 
the coherent fraction / — 1 — Nmc/N . The experimental 
values for coherent fraction were measured from interfer- 
ence fringe constrast. 

We take all parameters directly from the published val- 
ues for the experiments, with the exception of the gas 
temperatures. In the regime of very high condensate frac- 
tion, quantum gas thermometry is experimentally very 
difficult, and only a rough upper bound on the experi- 
mental temperature can be estimated {T/T^ < 0.36, for 
example, in [13|). We have therefore taken temperature 
as a fitting parameter, tuning it so that our theoretical 
curve hits the first experimental data point (i.e. the one 
for lowest lattice strength). This may be considered an 
application of 'coherence thermometry': we use ideal gas 
thermodynamics in the lattice to infer the gas temper- 
ature from interference contrast measurements. We as- 
sume that the gas temperature remains constant, within 
the same experimental data set, for all other barrier 
heights. As already mentioned, this isothermal assump- 
tion is rather crude, and could certainly be improved 
upon by considering detailed quantum kinetic studies of 
the effects of changing lattice strength. The fact that in- 
terpretation of the data in line with our isothermal the- 
ory curves is nevertheless surprisingly straightforward is 
a significant point, therefore, even if it should turn out 
to be co-incidental, because it bears on the sensitivity of 
this type of experiment to quantum gas kinetics. We will 
discuss this issue further in our final section. 

The various experiments we consider span the range 
of effective dimensionalities. In Fig. 4 is shown the case 
where the lattice heights Vi in all 3 directions are kept 
equal as they are ramped up, so that the gas is always 
effectively three-dimensional. The data for this case are 
taken from [l3|. Our fitted temperature is T = AbnK — 
0.29Tp, where Tq is the critical temperature for BEG 
in the (rather elongated) experimental trap, which had 
{a;^,o, ojyfi^^.A = 2vr X {20, 120, 120} Hz. 

In Fig. 4 we see that in the 3D geometry the noninter- 
acting, isothermal gas's coherent fraction changes very 
little for lattice heights from 2Er to AE^. For stronger 
lattices than this, it monotonically decreases, until about 
lO-E^; but at still higher lattice strengths, the coher- 
ence of the non-interacting gas actually increases slightly. 
This is because the stronger harmonic confinement of 
the gas is slightly over-compensating for the compres- 
sion of the density of states in the stronger lattice. As 
one would expect, the tight-binding approximation ap- 
proaches the full theory quite closely for strong lattices. 
The effective mass approximation becomes quite poor for 
stronger lattices, because it ignores the large band gap 
which is opening up. Since it includes thermal popula- 
tion of within-gap modes that do not actually exist, it 
greatly exaggerates the non-condensed population and 
hence underestimates the coherent fraction. 

The theoretical curve for the isothermal ideal gas has 
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FIG. 4: Coherent fraction of ideal Bose gas (solid line) as a 
function of optical lattice height Vi = V, with all three lattice 
components equal in strength. The dashed and dot-dashed 
lines show the results under the tight-binding and parabolic 
approximations described in the text. The dotted line de- 
scribes isentropic evolution. Points are experimental data 
taken from [13]. 



FIG. 5: Coherent fraction of the ideal Bose gas (solid line) as a 
function of transverse optical lattice height Vu-ans- The lattice 
strength along the long axis of the sample is held constant at 
Vax = 8Er. The dashed and dot-dashed lines give the same 
quantity under tight-binding and parabolic approximations, 
respectively. Points are experimental data taken from [32] . 



an excellent qualitative and quantitative agreement with 
experimental data for lattice heights up to ^ SE^, but 
thereafter the observed coherence of the real, interact- 
ing gas drops sharply. The fair agreement between the 
observations and the ideal gas in effective mass approx- 
imation is clearly co-incidental, since the drop in the ef- 
fective mass approximation curve at this point is an error 
in its representation of the ideal gas. As the lattice be- 
comes stronger, tunneling rates between lattice sites drop 
sharply, so that interactions become more and more im- 
portant in comparison with kinetic energies, while the 
atomic interactions also become gradually greater as the 
gas becomes more tightly confined within each site. The 
ideal gas approximation gradually worsens, then breaks 
down abruptly at the Mott transition. 

For this experimental case we have also computed the 
ideal gas coherent fraction assuming isentropic evolution, 
using a straightforward but computationally much more 
demanding version of our local density approximation. 
(Since entropy is an extensive variable, holding it con- 
stant is a non-local constraint.) As comparison between 
the solid and dotted curves of Fig. 4 shows, the dis- 
tinction between isothermal and isentropic evolution is 
not dramatic for lattice strengths below those at which 
the ideal gas model itself breaks down. In detail, how- 
ever, the isentropic dotted curve seems to do somewhat 
more poorly than the isothermal solid curve in fitting the 
'shoulder' or 'kink' feature of the data, which although 
it is a small effect, is common to all currently available 
experiments. (This 'coherency saturation kink' will be 
discussed further below.) This fact would seem to be 
worth noting for future non-equilibrium studies of dilute 
Bose gases in time-dependent lattices. 

In Fig. 5 we compare our ideal gas theory with exper- 



imental data reported in j33j . In this case an effectively 
two dimensional situation was realized, by ramping up 
one component of the lattice to the rather high value 
of 8Er, and subsequently raising the lattice components 
in the other two orthogonal directions together, keep- 
ing them equal. Here our fitting happens again to find 
T = AdnK ~ 0.29r^, while all the other parameters are 
as in [33]. Comparison shows that the ideal gas model 
still works very well for the same lattice height range 
of 2 to 8 i^r- The effective mass approximation is poor 
in this case because the lattice is always quite strong in 
one direction. A mixed approximation would presumably 
be better, assuming tight binding in one direction and a 
parabolic spectrum in the other two. 

A quasi-one-dimensional version of this problem has 
also been realized [3J], and is shown in Fig. 6. The 
trapped gas was in this case loaded into a strong two- 
dimensional lattice AOEr high, and a weaker lattice po- 
tential along the axial direction Vax was raised to various 
much lower heights. In this case we see excellent agree- 
ment between the isothermal ideal gas theory and the 
experimental data, albeit in a regime where the varia- 
tion is not dramatic in either one. To fit the initial ex- 
perimental point, however, we have had to take a quite 
high temperature, nearly twice what the initial critical 
temperature would be in the same harmonic trap with- 
out the lattice: T = 193nK = 1.94Tg. (This harmonic 
trap had {w^,o, i^y,o,^zM} = 27r x {8.7, 90, 90} Hz.) With 
the strong transverse lattice and initial axial lattice at 
Vax = 2Er, however, the ideal gas critical temperature 
for BEC is actually 250nK, so our fitted temperature was 
still well below the actual Tc- It is high enough, though, 
that the thermal occupation of the second axial band is 
non-negligible, and so in this case we have included it in 
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FIG. 6: Coherent fraction of the ideal Bose gas (solid line) 
as a function of axial optical lattice height Vax- The lattice 
height in the other two directions is kept fixed at Vtrans ~ 
40-Br. Points are experimental data taken from [33]. 



our full theory curve. 

The temperature of the initial condensate in this ex- 
periment, before loading it into the lattice, was reported 
to be only around 5QnK. It is unclear, therefore, whether 
heating during this loading process has really raised the 
temperature by a factor of four, or whether many-body 
effects are already suppressing coherence at weak axial 
lattices because of the tight transverse confinement, or 
whether the high gas compression of the transverse lattice 
has simply raised the rate of scattering among atoms dur- 
ing the time-of-flight imaging process enough to degrade 
the observed interference patterns. The simplest expla- 
nation of Fig. 6, which is that turning on the strong trans- 
verse lattice had indeed considerably heated the samples, 
is by no means extremely implausible, however. The ev- 
idence from coherence for quantum many-body effects, 
as opposed to ideal gas thermodynamics, would seem to 
be weaker in the one-dimensional geometry than in the 
higher dimensionality cases. 



A. Coherency saturation 

In all three cases discussed above, the initially falling 
theoretical curves level out, and then start to increase 
at some value of the potential height. This appar- 
ently counter-intuitive behavior of the non-interacting 
gas is in fact easily explained, as being due to the co- 
incidental effective strengthening of harmonic confine- 
ment, along with strengthening of the lattice barriers, 
when the focused red-detuned laser intensities are raised. 
A stronger lattice means a narrower energy band, hence a 
greater density of states and larger non-condensed popu- 
lation, hence a lower coherent fraction. But a tighter 
trap at fixed (critical) chemical potential holds fewer 
non-condensed particles, and hence favors higher coher- 
ence. Which of these competing effects dominates de- 
pends upon the lattice and trap strengths and shapes. It 



FIG. 7: Coherent fraction as a function of optical lattice 
height under experimental conditions of [SJ. The arrow 
shows the point where the energy spectrum change due to 
harmonic trapping starts to dominate the change due to the 
periodic potential. The inset is data taken from ^JIi a^nd 
gives a numerical derivative of the experimentally measured 
visibilty curve as a function of lattice height. 



also depends upon temperature. For temperatures rang- 
ing from T = 30 - eOnK (which is T ~ 0.1 - 0.2T^ 
for the elongated trap of (S^, with {ujxfi,ujyfi,uJzfi} ~ 
27rx {20, 200, 200} Hz), the lattice height at which the lo- 
cal coherency minimum occurs ranges over V ~ 14 — 8i?r. 
(Higher temperature brings the minimum to lower lattice 
heights.) Even for the ideal gas the quantitative details 
in this effect are involved, but the two competing factors 
are qualitatively simple and clear. 

It seems worth bearin g th em in mind when investi- 
gating observations, as in 35], of 'kinks' in the coherence 
fraction curve for the real gases N lattice heights at which 
the steady fall of coherent fraction seems to pause, be- 
fore resuming. The fact that the coherency does resume 
falling, instead of rising as in the ideal case, is evidently 
due to interparticle interactions. And the true explana- 
tion of these coherency 'kink' observations may indeed 
have to do with the shifting of supcrfluid and Mott in- 
sulator shells with lattice height in the presence of the 
harmonic potential. But as Fig. 7 shows, when the num- 
ber of atoms is 3.6 x 10^ and the temperature is tuned 
to r = 30nK, the ideal gas coherency minimum is at 
l4:Er, the same location seen experimentally. The fact 
that only mild re-tunings of the temperature are needed 
to make the ideal gas coherence curve show similar kinks 
at exactly the same lattice heights indicates that coher- 
ence enhancement because of rising effective harmonic 
confinement may also be playing a significant role at this 
point in the experimental system. It is beyond the scope 
of this paper, however, to assess the impact of this effect 
beyond the ideal gas regime. 



VI. CONCLUSION AND OUTLOOK 

In this paper we have shown that optical lattices 
formed from red-detuned gaussian profile laser beams 
contribute an effective harmonic confinement to a 
trapped quantum gas, and that this effect needs to be 
taken into account in computing the coherence fraction 
of an ideal Bose gas in such a lattice at a given sub-critical 
temperature. It is in principle quite possible to remove 
this effect by altering the magnetic trap while raising the 
laser intensity, so as to hold the total harmonic confine- 
ment constant while the lattice is raised. We have shown 
that doing so would significantly simplify the interpreta- 
tion of experimental results. 

The good agreement between experiment and ideal gas 
theory that we have found up to threshold lattice heights 
supports the interpretation that many-body effects are 
slight up to that point, but that true many-body effects 
such as number squeezing and even the Mott transition 
do then appear, because ideal gas physics does not seem 
able to account naturally for the rather sharp drop in co- 
herence observed at these points. It is important to note 
that the mildness of the variation with lattice height of 
the ideal gas's coherent fraction is largely due to cancella- 
tion between the lattice's band compression and the con- 
comitant increase in harmonic trapping strength, as we 
have just summarized. Without this coincidental com- 
pensation through effective harmonic confinement, our 
Fig. 3 showed clearly that band compression at constant 
temperature can sharply lower the coherence of an ideal 
Bose gas in a strong optical lattice. It is therefore impor- 
tant to consider the effects discussed in this paper when 
inferring many body physics from coherency data. 

A weak point in our theory as it stands is certainly 
the assumption that the temperature does not vary with 
lattice height. If the lattice is ramped on slowly enough 
for the cloud to remain in global equilibrium, and if the 
gas temperature is fixed by evaporation and heating pro- 
cesses that are insensitive to the lattice barrier height, 
then the isothermal assumption may well be valid. De- 
termining precisely where observations depart from equi- 
librium ideal gas theory may also require more careful 
consideration, however. It seems quite plausible that our 
isothermal equilibrium curves are to some extent mim- 
icking the more realistic combination of adiabatic cooling 
with non-equilibrium entropy generation. And perhaps, 
near the threshold lattice height, adiabatic cooling of the 
gas is able to compensate for some decoherence due to 
incipient number squeezing. If that were so, many body 
physics could actually be emerging somewhat before the 
experimental coherence fell below our isothermal theory 
curves. Or perhaps heating does increase before quan- 
tum many body effects really predominate, so that the 
first fall in coherence involves only ideal gas dynamics, 
and the onset of many-body dynamics comes only later. 
These issues will require further study, even just within 
ideal gas dynamics. Extension into weakly interacting 
dynamics, with mean field theory and Bogoliubov quasi- 



particles, is then an obviously desirable next step. 

Finally, crude though its application in this paper has 
been, our temperature fitting procedure has raised the 
prospect of applying coherence thermometry in experi- 
ments on cold bosons in optical lattices, in regimes where 
normal time-of-flight thermometry becomes imprecise. 
Better theoretical calibration, taking into account the 
non-trivial kinetics issues described above, may make this 
a precise technique for future experiments. 



VII. APPENDIX 

In this Appendix we derive the transcendental equa- 
tion from which the dispersion relation in the uniform 
lattice may be determined numerically. 

The one-dimensional time-independent Schrodinger 
equation for the biparabolic periodic potential may be 
written 



where 
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where ^{X) is a confluent hypergeometric function. 




and thus the wavefunction in this region can be written 
as 



^I = Ci(pi(0 +C2(^2(?), 



(A4) 



whith unknown constant coefRcients Ci and 02- 

Correspondingly the linearly independent solutions in 
the " barrier- type" region 37r/2 < ^ < 57r/2 are 



with unknown constant coefficients ci and C2. 

Using the Bloch theorem [36( we can write the wave- 
function in the region 57r/2 < ^ < 77r/2 as 
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Now implementing the continuity requirements of the 
wavefunctions and their derivatives at the boundary 
points 37r/2 and 57r/2, we get a homogeneous system 
for the unknown coefficients, and from the requirement 
of nontrivial solution we obtain the dispersion relation 
Eq.(7), with 



G(^§-,^^=iMO^'iiO + ^iiO^'iiO) (A9) 

X (^2(0^2(0 +^2(0¥'2(0)I?=V2: 



with prime denoting the first order derivative. 
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